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Erra.ta.: 

1.  Equation  (1)  should  be  changed  to 

fi.x)dz  =  k(p  nxi)  +  +  ^(-/(l-,  )  +  /(!,)  +  /(jn-2)  -  /(*„ 

2.  Equation  (2)  should  be  changed  to 

3^(/)  =  '■("E  /(li)  +  +  /*  E  o,if{x„., )  -  fix,)) 

i=l  j=l 

3.  The  weights  tabulated  in  Section  7.1  are  — and  not  as  they  should  be. 


A  group  of  quadrature  formulae  is  presented  applicable  to  both  non-singular  functions  and 
functions  with  end-point  singularities,  generalizing  the  classical  end-point  corrected  trape¬ 
zoidal  quadrature  rules.  We  present  an  algorithm  for  the  construction  of  very  high-order 
end-point  corrected  trapezoidal  rules,  taking  advantage  of  functional  information  outside 
the  interval  of  integration.  The  scheme  applies  not  only  to  non-singular  functions,  but  also 
for  a  wide  class  of  functions  with  monotonic  singularities.  Numerical  experiments  are  pre¬ 
sented  demonstrating  the  practical  usefulness  of  the  new  class  of  quadratures.  Tables  of 
quadrature  weights  are  included  for  singularities  of  the  form  s(x)  =  /opdi]),  s(i)  =  |i|’^  for 
a  variety  of  values  of  A. 
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1  Introduction 


The  trapezoidal  rule  is  known  to  be  an  easy  and  numerically  stable  means  for  numerical  in¬ 
tegration.  If  a  function  is  periodic  and  analytic  on  the  interval  of  integration,  the  trapezoidal 
rule  converges  exponentially  fast  (see,  for  example,  [7]).  However,  for  non-periodic  functions 
the  trapezoidal  rule  is  second  order  convergent,  and  end-point  corrections  are  often  used  to 
improve  the  convergence  rate.  A  standard  end-point  corrected  trapezoidal  rule  is  fourth  order 
convergent,  and  is  given  by  the  formula 

Hx)dx  =  £  /(li)  +  +  /(I,)  +  /(!„-,)  -  /(!„)), 

(1) 


where,  h  =  {b  -  a)/{n  -  1)  and  —  a  +  ih  for  i  =  0,  l,2,...,n  -  1  (see,  for  example  [1]). 


More  recently,  the  Euler-Maclaurin  formula  is  used  in  [4]  to  obtain  a  high-order  end-point 
corrected  trapezoidal  rule  of  the  form 


TM)  =  E  /(*,)  +  +  *  E  -  /(«>)). 


n-2 

E 

i=:l  “ 

where  a  =  (otijOj,  ...,am)  are  coefficients  such  that 


(2) 


TO)  -  l'f{x)dx\<^ 
Ja  n 


c 

m  ' 


(3) 


for  some  c  >  0. 


The  scheme  of  [4]  provides  satisfactory  quadratures  upto  order  12;  for  higher  orders,  the  coeffi¬ 
cients  Q  grow  rapidly,  rendering  the  scheme  useless.  In  this  paper  we  develop  a  different  class  of 
end-point  corrected  trapezoidal  rules,  whereby  the  growth  of  correction  weights  is  suppressed, 
enabling  the  construction  of  end-point  corrected  trapezoidal  rules  of  arbitrarily  high  order  for 
non-singular  functions. 

In  [5],  end-point  corrected  quadrature  formulae  are  developed  to  approximate  definite  inte¬ 
grals  of  singular  functions  /  :  [a,  6]  — ►  of  the  form 

f{x)  =  <l>(x)s(x)  +  (4) 

and 

f{x)  =  4>ix)s{x),  (5) 

where  a  <  0  <  b,  <f>{x),‘fl>(x)  €  c*[a,6],  and  s{x)  €  c[a,h]  is  an  integrable  function  with  a 
singularity  at  0.  The  procedure  developed  in  [5]  provides  satisfactory  quadratures  only  upto 
order  4;  for  higher  orders,  the  quadrature  weights  grow  rapidly,  also  rendering  the  scheme 
useless.  In  this  paper  we  construct  a  different  class  of  end-point  corrected  trapezoidal  rules. 
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whereby  the  growth  of  quadrature  weights  is  partially  suppressed  for  functions  of  the  form 
(4),  obtaining  useful  quadratures  of  order  upto  12;  and  completely  suppressed  for  functions 
of  the  form  (5),  providing  quadratures  of  arbitrarily  high  order.  Moreover,  we  obviate  the 
programming  inconvenience  associated  with  the  procedure  developed  in  [5],  which  requires  that 
functional  information  be  tabulated  on  a  grid  finer  than  that  required  for  the  uncorrected 
trapezoidal  rule. 

Remark  1.1  The  approach  of  this  paper  is  somewhat  related  to  that  of  [9].  However,  [9] 
constructs  quadratures  in  higher  dimensions,  and  these  quadratures  are  of  relatively  low  order. 
In  this  paper,  we  construct  one-dimesional  rules  of  very  high  order.  Furthermore,  most  rules 
of  this  paper  are  “standard”  in  the  sense  that  the  correction  coefficients  do  not  depend  on  the 
number  of  nodes  in  the  trapezoidal  rule  being  corrected,  or  on  the  sampling  interval. 

2  Mathematical  Preliminaries 

In  this  section  we  summarize  some  well-known  results  to  be  used  in  the  reminder  of  the  paper. 
Lemmas  2.1,  2.2  and  2.3  can  be  found,  for  example,  in  [1]. 

Definition  2.1  Suppose  that  a,b  are  a  pair  of  real  numbers  such  that  a  <  b,  and  that  n  >2  is 
an  integer.  For  a  function  f  ;  [a, 6]  —*  R^,  we  define  the  n-point  trapezoidal  rule  Tn{f)  by  the 
formula 

«=0  ^ 

with 

=  (6  -  a)/(n  -  1).  (7) 


The  following  lemma  provides  an  error  estimate  for  the  approximation  to  the  integral  given  by 
the  trapezoidal  rule.  Along  with  Lemma  2.2,  it  can  be  found,  for  example  in  [1]. 


Lemma  2.1  ( Euler- Maclaurin  formula)  Suppose  that  a,b  are  a  pair  of  real  numbers  such  that 
a  <  b,  and  that  m>  I  is  an  integer.  Further,  let  Bk  denote  the  Bernoulli  numbers 


B2 


1 

6’ 


(8) 


jy/ €  c^”*+2[a,6]  f  has  2m 2  continuous  derivatives  on  [a,b]),  then  there  exists  a  real 

number  with  a  <  ^  <b,  such  that 


f 


nx)dz  =  TM) + p  - 


f2m+2(c\ 
(2m  +  2)!  ^ 


(9) 
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The  following  well-known  lemma  provides  an  error  estimate  for  Lagrange  interpolation. 

Lemma  2.2  (Lagrange  interpolation  formula)  Suppose  that  a,  b  are  a  pair  of  real  numbers  such 
that  a  <  b,  m  >  3  be  an  odd  integer,  and  f  €  c”*[a,6].  Furthermore,  let  h  be  defined  in  (7),  and 
f  be  tabulated  at  equispaced  points,  xjt  =  kh.  Then  for  any  real  number  p  there  exists  a 

real  number  f ,  —mh  <  ^  <  mh,  such  that 


(m-l) 


fiXQ  +ph)=  A^(p)f(x/t)  +  Rm-l , 


^  2 


with 


and 


m-l 


m  —  1 


m— 1 


(10) 


(11) 


(12) 


Lemma  2.3  If  f  :  [a,6]  — ►  iZ*  is  a  function  satisfying  the  conditions  of  Lemma  2.2,  and  the 
coefficients  D!f),  are  given  by  the  formula 


=  ii5rii)(^r(P)) 

(13) 

then 

m— 1 

(14) 

for  any  m,i  such  that 

<k<  andl<k<  iSfi. 

Proof.  The  proof  is  as  an  immediate  consequence  of  (10)  and  (13). 

□ 

Lemma  2.4  Suppose  that  m,l,k  are  integers,  and  the  coefficients  a'j^i  are  defined  by  the  recur¬ 
rence  relation 


=  1, 

(15-a) 

®1,2 

=  1, 

(15-b) 

-  (*  -  *  i  -1- 

(15-c) 

-  “fc,/-2  (  2  1 

(15-d) 

3 


with  a^i  =  0,  for  all  k  <  0,  or  I  <  0,  or  m  <  1.  Then 


ATip)  = 


»  1  m  — « 

“5- 

zJ.  _  Jbv  ^  ’ 


/or  any  odd  m  >  3,  1  <  /:  <  and  v4][*(p)  w  defined  by  (11). 
Proof.  Due  to  (11), 


AD{p)  = 


(-!)■ 


'-+k 


where 


(V+A:)!(2!fi-ik)! 


cr(p), 

-0- 


Thus  it  is  sufficient  to  show  that 


m»-»l 

2 


cir(p)  =  E  “wp'- 

Isl 

This  will  be  shown  by  induction.  Indeed,  if  m  =  3  then,  due  to  (18), 

Cf{p)  =  p^+Py 


which  is  equivalent  to  (15-a),(15-b). 

Assume  now  that  for  some  m,k  such  that  — <  k  < 

m— I 

cr(p)=  Eo&p'- 

(=1 

Combining  (18)  and  (21),  we  have 

cr”(p)  =  (p+^)(p-^)i:»up' 
^  ^  1=1 
m— 1 

^  1=1 

m— I 

=  E»Sp'*’-(^)=E»&p'. 
/=!  ^  /=! 

which  is  equivalent  to  (15-d). 

Now,  assume  that  for  some  k 

cr‘(P)=E“S"v'- 

1=1 


(16) 


(17) 

(18) 


(19) 


(20) 


(21) 


(22) 
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(23) 


Combining  (23)  and  (18),  we  have 


Cr+"(P)  =  (P-*)(P-(A:+I))53a*f+V 

/=i 

=  (pH  P  -  +  *))  52  V' 

/=i 

l=l  /=!  /=1 


which  is  equivalent  to  (15-c). 


Lemma  2.5  Suppose  that  Tn>  3  is  odd.  Then, 

=  (2i^  +  A)!(2i^-;b)!“^"‘-'^^*  " 

for  any  k,  i  such  that  <  k  <  and  I  <i  < 

with  the  coefficients  d^i  defined  by  the  recurrence  relation  in  Lemma  2.4. 
Proof.  Substituting  (16)  into  (13),  we  immediately  obtain 

nm  _  K~^)  ^ _ ^  m  J  |p=0 

’■*  +  k)\i^  -  it)!  9p(2*-J)  ^  ' 

(21fi  +  jfc)!(2ifi  -  ~ 


The  following  six  lemmas  provide  identities  which  are  used  in  the  proof  of  Theorem  3.1. 
Lemma  2.0  If  k  >2  is  an  integer  and  a^j  is  defined  in  Lemma  2.4- 

<  ia+2)!-45«  I,  (27) 

for  all  I  =  1,2,  ...,2fc-  3. 

Proof. 

If  A:  =  2,  and  /  =  1  then  |  (1)!  •  j  |=  2,  |  (3)!  •  02,3  i=  12,  and  therefore  (27)  is  obviously  true. 
Now,  assume  that 

KOI-^fA'I  <  I  ('  +  2)!- -St’,  I,  (28) 

for  some  A:  >  2  and  all  /=  1,2,  ...,2A:  -  3. 

Now,  due  to  (15-a),  (15-b),  (15-c),  and  (15-d), 

(')!  •  4+u  =  (((*+!)-(*+  l)’)og«  +  a|‘i>  +  alf_\)  ■  (I!),  (29) 


5 


and 


('  +  2)!  ■  »“15+2  =  m  +!)-(*  +  +  agS  +  »”*')  ■  ('  +  2)!- 

Finally,  combining  (28),  (29),  and  (30)  we  easily  obtain 

<  l('  +  2)!-a2Sj+2l. 

for  all  1  =  1,2,  ...,2k  -  1. 


(30) 

(31) 
□ 


Lemma  2.7  If  k>  2  is  an  integer,  and  is  defined  in  Lemma  2.4  then 


I  (/)!•<,  I  <  l(/  +  2)!-</+2l. 


(32) 


for  all  m  >  2k  1  and  1=1,2, ...,  21?  -  3. 

Proof.  Lemma  (2.6)  establishes  the  base  case,  i.e.,  that  (32)  is  true  when  m  =  2fc  +  1. 


Now,  assume  that 

I  (/)■'•</ I  <  l(/  +  2)!-a]^,+2l,  (33) 

for  some  odd  m  >  2A:  +  1,  and  all  1  =  1,2,  ...,2A:  —  3. 

Now,  due  to  (15-a),  (15-b),  (15-c),  and  (15-d), 

(()!  •  =  (aj;,.;  -  (^)’aa)  •  ((!),  (34) 

and 

('  +  2)!  •  =  (“&  -  •  ('  +  2)!-  (35) 

Finally,  combining  (33),  (34),  and  (35)  we  easily  obtain 

(0!-|<r'l  <  l(^  +  2)!-a;jVu+2l’  (36) 

for  all  1  =  1,2,. ..,21:- 3.  □ 

Lemma  2.8  If  m,k  are  integers  such  that  m>3  is  odd,  and  ^  k  <  then 

I  (1)!  •  aT,^  |<|  (3)! .  aT,,  |<|  (5)!  •  a^,  ]<  ...  <1  (m  -  2)!  •  1  .  (37) 

Proof.  This  Lemma  follows  directly  from  Lemma  2.6  and  Lemma  2.7  □ 


Lemma  2.0  If  m>Z  is  odd,  then 


(m-  l)(m-2)!  ^  (27r)'"-* 
2((=^)!)2  4 


(38) 


6 


I 


Proof.  If  m  =  3  then  obviously  1  < 
Now,  assume  that  for  some  odd  m  >  3, 


(m  — l)(m-2)!  ^  (2^)”*  ' 

2((2^)!)2  <  — r~ 


Obviously, 

(m-b  l)(m) 

(il2^)2 

and  combining  (39)  and  (40)  we  obtain 


Am 

(m  +  1) 


<  (25r)^, 


(m+  l)(m)(m-  l)(m-2)!  (27r)*"  *.0  x2 

(("»+>) )2  2((2^)!)2  4  ^  ’ 

which  is  equivalent  to 

(m  +  l)(m)!  ^  (2x)’"+^ 

2((Ma)!)a  <  4  ■ 

Now,  the  conclusion  of  the  lemma  is  an  immediate  consequence  of  (39)  and  (42). 


(39) 


(40) 


(41) 


(42) 

□ 


Lemma  2.10  //m  >  3  is  odd  then 


1  I  < 


(2ff) 


,m-l 


m-il 


for  any  fc,  i  such  that  <  k  <  and  1  <  i  <  — y-* 

Proof.  Combining  Lemmas  2.5,  2.6,  2.7,  and  2.8,  it  is  easy  to  see  that 

1  1  <  1  I  <  -  <  1  I  • 

Consequently,  it  is  sufficient  to  show  that 

(27r)”‘-^ 


Dm 


m  — 1  r 
7 


l< 


First  we  observe  that  (obviously)  for  any  k  such  that  - 


k{m  -  2)! 


2  ^  ^  2  ■ 


_ (m-  l)(m-2)! 

(=^  +  k)\{^--^  -  k)\  2((2ifl)!)2  ■ 

Then,  we  combine  (15-a),  (15-b),  (15-c),  (15-d),  and  (25)  to  obtain 

1  T-.m  I  ~  2)! 

I  1=  (a^  +  ifc)!(iSfi-ifc)!- 

Now,  (45)  follows  immediately  from  the  combination  of  (46),  (47),  and  Lemma  2.9. 


(43) 


(44) 

(45) 

(46) 

(47) 

□ 


7 


Lemma  2.11  For  any  I  >  I  the  Bernoulli  number  B21  satisfies  the  inequality 


I  -V  I  <  — i _ 

'  (2/)!  '  (2t)2'* 

Proof.  As  is  well  known  (see  for  example,  [1]),  for  any  /  >  1 

(-iy-»2(2/).«^  1 

(27r)2' 

and 

00  . 

H  p7 

fc=i  * 

Now,  the  conclusion  of  the  lemma  is  an  immediate  consequence  of  (49)  and  (50). 
The  proof  of  the  following  lemma  can  be  found  in  [5]. 


(48) 


(49) 


(50) 

□ 


Lemma  2.12  Suppose  that  m  >  1,  s  €  c”*(0, 1)  possesses  a  finite  integral  on  the  interval  [0, 1], 
and  that  is  monotonic  in  some  neighborhood  ofQ.  Then  the  product  x-s(i)  is  bounded  on 

[0, 1).  Suppose  further  that  w  6  c'"[0, 1]  is  such  that  u;(0)  =  u;’(0)  =  u)”(0)  =  ...  =  u)^"*^(0)  =  0. 
Then  the  function  0(i)  =  s(x)  •u;(z)  is  defined  on  the  closed  interval  [0, 1],  and  V'(O)  =  ^'(0)  = 
V»"(0)  =  ...  =  tl’^"‘)(0)  =  0. 


3  End-point  Corrections  for  Non-singular  Functions 

3.1  End-point  corrected  trapezoidal  rules 

While  the  authors  have  failed  to  find  the  contents  of  this  section  in  the  literature,  it  is  an 
immediate  consequence  of  well-known  facts  from  classical  analysis.  We  present  it  here  for  com¬ 
pleteness,  and  because  we  found  the  resulting  high-order  quadrature  rules  quite  useful  (see 
Section  7.1). 


Suppose  that  n,  m,  are  a  pair  of  integers  with  m  >  3  and  odd,  and  n  >  2.  Further,  suppose  that 
a,b  are  a  pair  of  real  numbers  such  that  a  <  b,h  =  (6-a)/(n-l),  and  /  :  [a  —  mh,b+mh]  -* 
is  an  integrable  function.  We  define  the  corrected  trapezoidal  rule  Tpm  for  non-singular  func¬ 
tions  by  the  formula 

TSAf)  =  Tnif)  +  h  {fib-¥kh)-f{a-\-kh)W.  (51) 

- - j 

The  real  coefficients  are  given  by  the  formula 


(=i 


(21)!  ' 


(52) 
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where  are  defined  in  (13)  (also,  see  (25))  and  B21  are  the  Bernoulli  numbers. 

We  will  say  that  the  rule  Tpm  is  of  order  m  if  for  any  /  6  £"*[0  -  mh,b  +  mh],  there  exists  a 
real  number  c  >  0  such  that 


lr^"m(/)-  tfix)dx\<^.  (53) 

Tl 

Theorem  3.1  If  m>  Z  is  an  odd  integer  then  for  any  k  such  that  — <  k  < 

(54) 

where  the  coefficients  are  defined  in  (52). 


Proof.  Combining  Lemma  2.10  and  Lemma  2.11  we  immediately  observe  that 


and  hence 


DJX 

(2/) 


1, 


m— 1 


m  -  1 
2 


(55) 


(56) 

O 


Remark  3.1  A  somewhat  more  involved  argument  shows  that  in  fact  <  1  for  all  k,m\ 
empirically  this  can  also  be  seen  from  the  tables  in  Section  7.1  below.  However,  for  the  purposes 
of  this  paper  (56)  is  sufficient. 


Theorem  3.2  Suppose  that  m,n  are  a  pair  of  integers  with  m  >  3  and  odd,  and  n  >  2. 
Further,  suppose  that  a,  b  are  a  pair  of  real  numbers  such  that  a  <  b.  Then,  the  end-point 
corrected  trapezoidal  rule  T^m  is  of  order  m,  i.e.,  for  any  f  :  [a  -  m/i,6  +  mh]  —*■  such  that 
f[a  -  mh,h  +  mh]  €  c’”[o  -  mh,b-\-  mh],  there  exists  a  real  number  c  >  0  such  that 

\T^M)-  t  f(x)dx\<^.  (57) 

Ja  ^ 

Proof.  Combining  (52)  and  (51),  we  obtain 


TSM) 


m— 1  m—-l 

— ^  ^  B 

TM)  +  h  £  (/(6+<:A)-/(a  +  *4))£-^ 

/=!  ^ 

K-  2 

^  ^  +  kh)  -  f{a  +  kh)) 

InU)  +  ^  (21)! 


(58) 
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Combining  (14)  and  (58),  ve  have 


m-»l 


nAD  =  r„(/)  +  E 

Finally,  combining  (59)  with  Lemma  2.1,  we  observe  that  for  some  a  <  ^  <  6, 


TMf)=  I  /(x)dx  +  2i2^!.-^>  +  ^^r(0, 

j(X  7Tl» 


and  the  theorem  immediately  follows  from  (60). 


(59) 


(60) 

□ 


Remark  3.2  It  is  easy  to  see  that  for  m  >  3  and  odd,  and  any  k  such  that  — <  k  < 

^ro  ~  ®  (13)),  and  hence  and  =  0  (due  to  52). 

Now,  instead  of  (51)  one  could  define  the  end-point  corrected  trapezoidal  rule  by  the  formula 

m— 1 

T^Af)  =  Tnif)  +  hfi  (Kb  +  kh)  -  f{b  -  kh)  -  /{a  +  kh)  +  /(a  -  kh))p)^.  (61) 

k=l 

4  End-point  Corrections  for  Singular  Functions 

In  this  section  we  construct  a  group  of  quadrature  formulae  for  end-point  singular  functions, 
generalizing  the  classical  end-point  corrected  trapezoidal  rules.  The  actual  values  of  end-point 
corrections  are  obtained  for  each  singularity  as  a  solution  of  a  system  of  linear  algebraic  equa¬ 
tions.  All  the  rules  developed  in  this  section  are  simple  extensions  of  the  corrected  trapezoidal 
rule  Tpm  developed  in  the  preceding  section. 

A  right-end  corrected  trapezoidal  rule  T^pm  is  defined  by  the  formula 

rV(/)  =  /»(^^^^  +  E/(xi))  +  hf2{fib+kh)-f{b-kh)W,  (62) 

^  1=1  fe=i 

where  /(O,  b  -f-  mh]  -*  is  an  integrable  function,  n,  m  are  a  pair  of  natural  numbers  with 
m  >  3  and  odd,  the  coefficients  /JJJ*  are  given  by  (52),  and 


Xi  =  ih.  (63) 

We  will  say  that  the  rule  is  of  right-end  order  m  >  3  if  for  any  /  €  m/i]  such 

that  /(0)[214z  =  /'(O)  =  ...  =  /^”*^(0)  =  0,  there  exists  c  >  0  such  that 

|rA!-(/)-//Wi<»l<^.  (61) 
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It  easily  follows  from  Theorem  3.2  that  is  of  right-end  order  m. 
Similarly,  a  left-end  corrected  trapezoidal  rule  is  defined  by  the  formula 


+  +  h^(-/(-b  +  kh)  +  /(-b-kh))j3r,  (65) 


where  /[—b  —  mA,  0)  — >  is  an  integrable  function,  n,  m  are  a  pair  of  natural  numbers  with 

m  >  3  and  odd,  the  coefficients  /3J^  are  given  by  (52),  and  /i,Xi  are  defined  by  (63).  We  will 
say  that  the  rule  is  of  left-end  order  m  >  3  if  for  any  /  €  6  —  m/i,0)  such  that 

/(O)  =  /'(O)  =  ...  =  /^’"^(O)  =  0,  there  exists  c>  0  such  that 


It  also  easily  follows  from  Theorem  3.2  that  T£jgm  is  of  left-end  order  m. 

Suppose  now  that  the  function  f{~kh,b  +  mh]  — ►  is  of  the  form 

f{x)  =  4>ix)s{x)  -f  V»(x),  (67) 

with  6  c*'(-Arh,5-f-mh],  and  s  G  c{-kk,b'\-mh]  an  integrable  function  with  a  singularity  at 
0.  For  a  finite  sequence  a  =  (a_jt,a_(fc_i),a_i,ai,...,a^)  and  T^m  defined  in  (62),  we  define 
the  end-point  corrected  rule  by  the  formula 


r.V(/)  =  r^«(/)  +  h  E  a,/(x,),  (68) 

j=i~k,j^O 

with  h  =  6/(n  -  1),  Xj  =■  jh. 

We  will  use  the  expression  T^pm  with  appropriately  chosen  a  as  quadrature  formulae  for 
functions  of  the  form  (67),  and  the  following  construction  provides  a  tool  for  finding  a  once 
pm  _  is  given,  so  that  the  rule  is  of  order  k,  i.e.,  there  exists  a  c  >  0  such 

that 

\T£pmif)-J^f{x)dx\<-^.  (69) 

For  a  pair  of  natural  numbers  k,  m,  with  fc  >  1  and  m  >  3  and  odd,  we  will  consider  the  following 
system  of  linear  algebraic  equations  with  respect  to  the  unknowns  a",  with  j  =  0,  ±1,  ±2, ...,  ±k: 

_ _ t  .--in  ”  JO 


j=-k,i^0 


for  i  =  1,2,  ...,fc,  and 


X;  =  i  J‘x<-i-'s{x)dx  -  »(!)), 


}=-k,j^0 


11 


for  i  =  fc  +  l,fc  +  2,...,2ky  with  h  =  6/(n  —  1),  Xj  =  jh  and  T^pm  defined  by  (62).  We  denote 
the  matrix  of  the  system  (70),  (71)  by  ^4"*,  its  right-hand  side  by  and  its  solution  by 
On  =  (o"fc>o”(fc_i)>-">olli,Q!i,...,0!j).  The  use  of  expressions  T^npm  as  quadrature  formulae 
for  functions  of  the  form  (67)  is  based  on  the  following  theorem. 

Theorem  4.1  Suppose  that  a  function  s  :  {—kh,b  +  mh]  -+  is  such  that  s  €  c'’{—kh,b  +  mh] 
and  s*  is  monotonic  on  either  side  of  0.  Suppose  further  that  the  systems  (70),  (71)  have 
solutions  (Q”j^,a”j^_jj,Q”i,Oj,...,a)J)  for  all  sufficiently  large  n,  and  that  the  sums 

£  (aj)=  (72) 

i— *jilO 

are  bounded  uniformly  with  respect  to  n.  Finally,  suppose  that  the  function  f  :  (—kh,b  +  mh]  -+ 
is  defined  by  (67).  Then,  there  exists  a  real  c  >  0  such  that 

IT^npmif)-  f'fix)dx\<^  (73) 

Ja  ^ 


for  all  sufficiently  large  n. 

Proof.  Applying  the  Taylor  expansion  to  the  function  /  at  x  =  0  we  obtain 

fix)  =  P(/)(x)  +  Rki<k)(x)six)  +  R/tW, 

t=0  '  '  «=0 

and  Rk{^),  Rki^)  are  such  functions  [-kh,b+  mh]  — ►  R^  that 

si(«(o)  =  Riwm = ...  =  4‘’(«(i))  =  0. 

R'imo) = ^(V'Ko) = ... = RfHm) = 0- 

Substituting  (74)  into  (73),  we  obtain 

I  -  f  mdx  I  <  I  -  [  i>(/)(x)rfi  I  + 

Jo  JO 

I  TSx,0j{Rk{4>)  •  S)  +  RkW)  -  A(iifc(^(x))s(x))(x)  -h  Rki^{x)))dx  I . 

JO 

Due  to  (70),  (71) 

TSnpJPU))  -  [  P(/)(x)dx  =  0, 

Jo 


(74) 

(75) 

(76) 

(77) 

(78) 

(79) 
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and  we  have 


\TSn0jf)-  f'fix)dx\< 

Jo 

I  RT^^,{s .  Rkm  -  .  Rk{m^)dx  1 

+  I  RT^^.iRkW)  -  f^(Rk(m^)dx  I 
+  I  i:]UM<l>)ms{jh)a])  +  {Rkmjh)a])  I .  (80) 

Due  to  (77)  and  (64),  there  exists  Ci  >  0  such  that 

I  RT^^kiRkW)  -  j\Rkim=^)dx  |<  (81) 

Combining  (76),  (64),  and  Lemma  2.12  we  conclude  that  for  some  C2  >  0 

I  RT^^,(s  •  R,(d>))  -  j\s  .  R,{m^)dx  |<  (82) 

Finally,  combining  (76),  (77)  and  Lemma  2.12  we  conclude  that  for  some  C3  >  0, 

I  E  (i!j(«0)i)s(jfc)oJ)  +  (Jij(V.)(ift)aJ)|<g.  (83) 

Now,  the  conclusion  of  the  theorem  follows  from  the  combination  of  (81),  (82),  and  (83).  □ 


4.1  Convergence  Rates  for  Singularites  of  the  forms  |x|^  and  log[\x\) 

For  the  reminder  of  the  paper,  will  denote  functions  {-kh,b  +  Tnh]  — ►  R^  defined 

by  the  formulae 

<i>i{x)  =  x'-\  (84) 

for  i  =  l,2,...,k,  and 

4>i{x)  =  x'-'‘~^s{x),  (85) 

for  t  =  A:  +  1,  +  2, 21:.  The  following  lemma  is  a  particular  case  of  a  well-known  general  fact 
proven,  for  example,  in  [8]. 

Lemma  4.2  If  5(x)  =  x^  with  A  a  real  number  such  that  0  <  |A|  <  1,  then  the  functions 
<t>\i4>2i-—,d>2k  constitute  a  Chebyshev  system  on  the  interval (-kh,b  +  mh]  (i.e.,  the  determinant 
of  the  2k  X  2Ar  matrix  Bij  defined  by  the  formula  Bij  =  4>i{tj)  is  non-zero  for  any  2k  distinct 
points  on  the  interval  {-kh,b  +  mh]). 

Theorem  4.3  If  s{x)  =  with  0  <  jAj  <  1,  then  the  convergence  rate  of  the  quadrature  rule 
T^n0^  is  at  least  k. 
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Proof.  It  immediately  follows  from  Lemma  4.2  that  the  matrix  of  the  system  (70),  (71)  is 
non-singular.  We  rescale  the  system  (70),  (71)  by  multiplying  its  ith  equation  by  for 
t  =  1,2,  ...,fc,  and  by  ,  for  t  =  /b  -I- 1,«;  -I-  2,  ...,2k,  obtaining  the  system  of  equations 

T,  =  n(  t  (86) 

t  =  1,2,  ...,k,  and 

2  ^  -  r5^n.(x‘-*-^+")),  (87) 

j=-kj^0 

for  t  =  fc  +  1,  A  -I-  2,  ...,2k. 

We  will  denote  the  matrix  of  the  system  (86),  (87)  by  Bk,  and  its  right  hand  side  by  ZJ.  Ob¬ 
viously,  Bk  is  independent  of  n,  and  using  Theorem  3.2  we  observe  that  if  m  >  fc  ihen  \Z^\  is 
bounded  uniformly  with  respect  to  n.  Now,  due  to  Theorem  4.1,  the  convergence  rate  of  T^npm 
is  at  least  k.  □ 

The  proof  of  Theorem  4.3  can  be  repeated  almost  verbatim  with  s(x)  =  /o5(|x|),  instead 
of  s(x)  =  |x|^,  resulting  in  the  following  theorem. 

Theorem  4.4  If  s{x)  =  /o^di))  then  the  convergence  rate  of  the  quadrature  rule  is  at 

least  k. 

4.2  Asymptotic  behaviour  of  correction  coefficients  as  n  — »  oo 

An  obvious  drawback  of  the  expressions  T^pm  as  practical  quadrature  rules  is  the  fact  that 
the  weights  a"  =  (q"^,  ...,a"i,a",  ...,0^)  have  to  be  determined  for  each  value  of  n  by 
solving  a  system  of  linear  algebraic  equations.  For  singularities  of  the  form  s(x)  =  /oj(|x|), 
s(x)  =  Ixl'^  we  eliminate  this  problem  by  constructing  a  new  set  of  quadrature  weights  7  = 
(7_ife,7_(fc_i),  ...,7_i,7i,  -.,7*),  independent  of  n,  and  such  that  the  quadrature  rules  TJ^kpm  are 
still  of  order  not  less  than  k. 

Lemma  4.5  Suppose  that  P  =  is  such  that  the  right-hand  order  of  the 

quadrature  formula  TJ^m  is  m.  Further,  let  z  >  0  be  some  real  number.  Then  for  any  integers 
p,  q  such  that  p  <  q, 

I  -  [  *'«'*)  -  -  [  »'*)  1=  ).  (88) 

where  hp  =  b/{p  -  1),  and  hq  =  b/{q  -  1). 

Proof.  Due  to  Theorem  3.2,  there  exist  real  C|,C2  >  0  such  that 

(r^m(x')  -  /'**dx)  =  -hpj2  (89) 

•'0  j=-fc 
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and 


(90) 


-  f'x^dx)  =  C2h’^  -h,Yl  ilKY- 
Jo  j=-k 

Now,  combining  (89),  (90)  we  obtain 

I  ;AT(rJ5m(*')  -  /o -  /o‘ »''<<*)  I 


«p  j=-A:  is=-Jfc 

=  (ciV'"'  -  E  0)')-  (':*?■'■'  -  E  (■'■)’) 

j--k  3—-^ 

=  0(h^-^-^). 


Theorem  4.6  Suppose  that  k,  m  are  two  natural  numbers  such  that  k  <  m  —  1 
P  =  (/3{",/?J‘,,..,/3”_,J  is  such  that  the  right-hand  order  of  the  quadrature  T/j/jm  is  m. 

3 

further  that  s(i)  =  jx]^  with  0  <1  A  |<  1,  and  that  the  coefficients  (of!lfc,Q!!l(fc_i),a"i,c» 
are  the  solutions  of  the  system  (10),  (71).  Then 


I )  There  exists  a  limit 


limn-^ooOi  =  7i, 


for  each  i  =  1,2,  ...,2A:. 

2)  For  all  i  =  1,2,  ...,2fc, 

1  <  -  1=  0(^)- 

S)  7t  do  not  depend  on  m,  as  long  as  m  >  k  +  1. 

4)  The  quadrature  formulae  T!)pm  are  of  order  at  least  k. 


Proof.  Suppose  that  p,  q  are  two  natural  numbers,  and  p  <  q.  Obviously, 


a’  =  (Bi)-'ZJ, 

o<  =  (fl.r'ZJ, 
a” -a’  =  (B,)-‘(ZJ-Z|). 


Due  to  Lemma  4.5,  there  exists  c  >  0  such  that 


||ZJ-ZJ||< 


c 


(91) 
□ 

and  that 
Suppose 

? . «2) 

(92) 

(93) 

(94) 

(95) 
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and  by  combining  (94),  (95),  we  see  that  for  some  d  >  0 


Since  the  weights  a”  constitute  a  Cauchy  sequence,  they  converge  to  some  limit 
7  =  (7_*, 7-(ife-i).  7-i>  7i,  "M 7*),  which  proves  1,  and  2,  3,  4  follow  easily.  □ 

The  proof  of  the  following  theorem  is  a  repetition,  almost  verbatim,  of  the  proofs  of  the  Lemma 
4.5  and  Theorem  4.6 

Theorem  4.7  1/  under  the  conditions  of  Theorem  4-6  we  replace  s(x)  =  |x|^  with  s(x)  = 
/o^(|x|),  conclusions  1-4  remain  correct. 

For  singularities  of  the  form  li|^  and  /o^(|x|).  Theorem  4.6  and  4.7  reduce  the  quadratures 
to  the  more  “conventional”  form 


npr.if)-i-h  ^  7i/(®i).  (97) 


Remark  4.1  The  whole  theory  in  sections  4.1-4.2  has  been  constructed  for  functions  with  a 
singularity  at  the  left  end  of  the  interval.  Obviously,  an  identical  theory  holds  for  functions 
with  a  singularity  at  the  right  end  of  the  interval.  However,  in  all  formulae  the  expression  T^m 
has  to  be  replaced  with  (see  (62),  (65)). 


4.3  Central  Corrections  for  Singular  Rinctions 

In  this  section,  we  will  be  considering  functions  f[—b  —  mh,  0)  U  (0, 5  +  mh]  — ►  of  the  form 

/(x)  =  ^(x)s(x)  +  0(x),  (98) 

with  €  c'[-6  -  mh,  b  +  mh],  and  s  €  c[-6  -  mh,  0)  U  (0, 6  +  mh]  an  integrable  function  with 
a  singularity  at  0.  We  will  define  the  central-point  corrected  trapezoidal  rule 

l 

TJinp^if)  =  +  Tl0mif)  +  hJ2f^lifi^i)  +  /(a^-i)).  (99) 

3=1 

with  h,  Xj  defined  by(63),  defined  by  (52),  defined  by  (62)  and  (65)  respectively, 

and  /i"  =  an  arbitrary  sequence  of  length  i. 

We  will  use  the  expression  TJfnpm  with  appropriately  chosen  /i"  as  quadrature  formulae  for 
functions  of  the  form  (98),  and  the  following  construction  provides  a  tool  for  finding  /i"  once 
is  given,  so  that  the  rule  is  of  order  21,  i.e.,  there  exists  some  c  >  0  such  that 

I  i;.s"(/)  -  l<  -^r  (100) 
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For  a  pair  of  natural  numbers  /,m,  we  will  consider  the  following  system  of  linear  algebraic 
equations  with  respect  to  the  unknowns  /i”: 

^  (101) 

i=i  •'**’ 

for  i  =  1,2, and 


l 

E 

3=1 


.2«-2- 


=  t  x^'-^-^‘s(x)dx  -  Tk0„,{x^'-^-^‘s{x))  -  TEp^{x^'-^-^‘s{x)), 

J  —  h 


(102) 


for  i  =  /  +  1,  /  +  2, ...,  2/,  with  h  =  6/(n  —  1),  xj  =  jh. 

The  proofs  of  Theorem  4.8,  4.9,  and  4.10  are  almost  identical  to  those  of  Theorems  4.1,  4.3, 
and  4.6  respectively,  and  are  thus  stated  below  without  proof. 


Theorem  4.8  Suppose  that  a  function  s  :  [-h  —  mh,0)  U  (0,  i»  +  mh]  — >  B}  is  such  that  s  £ 
c^[-i>  —  mh,  0)  U  (0,  h  +  mh]  and  s*  is  monotonic  on  either  side  of  0.  Suppose  further  that  the 
systems  (101),  (102)  have  solutions  ,  •••»^")  for  all  sufficiently  large  n,  and 

that  the  sums 

E  w)'  ('“) 

i=-hj¥o 

are  bounded  uniformly  with  respect  to  n.  Finally,  suppose  that  the  function  f  :  [—6  —  mh,  0)  U 
(0,6+  mh]  -♦  B}  is  defined  by  (98).  Then,  there  exists  such  c  >  0  that 

|T;n^m(/)-  f' fix)dx\<^^  (104) 

Ja  ^ 

for  all  sufficiently  large  n. 


Theorem  4.9  If  s{x)  =  (ij’^  with  0  <  |A|  <  1,  ors(x)  =  /op(|x|),  then  the  convergence  rate  of 
the  quadrature  rule  TJ^np^  is  at  least  21. 

Theorem  4.10  Suppose  that  k,m  are  two  natural  numbers  such  that  k  <  m  —  1  and  that 
13  =  )  is  such  that  the  right-end  order  of  the  quadrature  Tupm  is  m,  and  the 

left-end  order  of  the  quadrature  Ti^m  is  m.  Suppose  further  that  s(x)  =  ji]^,  0  <]  A  |<  1,  or 
s{x)  =  lop(li|),  and  that  the  coefficients  » — *1^2)  ore  the  solutions  of  the 

system  (70), (71). Then 


1)  There  exists  a  limit 


=  fit, 


for  each  i  =  1,2,  ...,2h. 

2)  For  all  i  =  1,2,  ...,2k, 

ur-«i=0(;^)- 

S)  Pi  do  not  depend  on  m,  as  long  as  m  >  1  +  1. 

4)  The  quadrature  formulae  TJ^pm  are  of  order  at  least  21. 


(105) 


(106) 
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For  singularities  of  the  form  |x|'^  and  /o5(|i(),  the  Theorem  4.10  reduces  the  quadrature  to  the 
more  “conventional"  form 


f{x)dx « 

j-h 


4.4  Central  Corrections  for  Singular  Functions  f{x)  =  <^(x)s(i) 

In  this  section  we  construct  a  quadrature  formula  specifically  for  the  purpose  of  approximating 
definite  integrals  of  functions  of  the  form 

/(x)  =  ^(xMx),  (108) 

where  <^(x)  ;  c'’[-6  -  m/i,6  +  mh]  — ►  and  s  €  c[-6  -  mh,0)  U  (0,6  +  mh]  an  integrable 

function  with  a  singularity  at  0.  For  a  finite  sequence  p  =  (pQ,pi,P2t—>Pp)f  Tjipm,  T^pm 
defined  in  (62)  and  (65)  respectively,  we  define  the  corrected  trapezoidal  rule  T”„^n.  by  the 
formula 

T;^0mif)  =  n0M)  +  Tl0M)  +  h'£pli<i>ijh)  +  <i>(.-jh))).  (109) 

3=0 

For  integers  n,Tn,p  where  n>2,  p>l,m>3  and  odd,  we  will  consider  the  following  system 
of  equations  with  respect  to  the  unknowns  p"  =  (Po»Pi»P2»”M/»p): 

E  =  I  r  (®''‘'^(*))dx  -  r^n.(x2'-25(x))  -  r£pm(x2‘-2s(x)),  (110) 

where,  h  =  bf{n-  1),  xj  =  j/i,  and  t  =  1,2,  ...,p+  1. 

The  proof  of  the  following  theorem  is  almost  identical  to  the  proof  of  Theorem  4.1. 

Theorem  4.11  Suppose  that  n  >  2  is  an  integer,  and  h,Xi  are  defined  by  (6S).  Further, 
suppose  that  /(i)  =  4>{x)s{x)  where  4>  •  [-6  —  mh,b  +  mh]  — ►  R^,  and  s  €  c[— 6  -  mh,0)  U 
(0, 6  +  mh]  is  an  integrable  function  with  a  singularity  at  0.  Finally,  suppose  that  the  system 
of  equations  (110)  has  a  solution  (po>Pi» —»Pp)  fo^  ony  sufficiently  large  n  and  that  the  sums 
J2j=:o(Pj)^)  bounded  uniformly  with  respect  to  n.  Then  there  exists  a  real  c  >  0  such  that 

I  !?,(/) -/_/(*)■'*  l<  (“>) 

The  proof  of  the  following  theorem  is  almost  identical  to  that  of  Theorem  4.6,  and  is  omitted. 

Theorem  4.12  Suppose  that  s{x)  =  /op(|x|).  Then  for  all  n  >  2p,  the  system  (110)  has  a 
solution  p”^  =  (Po,Pi,p2>-.Pp)» 

Po  =  j  J  Jog{\x\)dx-T^mlog{\x\)-Tl0„log{\x\)-J2pj.  (112) 
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Pui  Ihermore,  there  exist  such  real  numbers  pi,P2i—iPp  ond  a  real  d>  0  such  that 


limn—roopj  —  pj^ 

and 

\p]-pj\<d-h^'^ 

for  all  j  =  1,2,  Finally,  there  exists  a  real  cq  such  that 

I  Po  -  +  O.blog{h)  -^Pj\<d‘  /i”*"'’ 

j=i 


(113) 

(114) 


(115) 


for  all  n  >  2p. 

Remark  4.2  Formulae  (114),  (115)  indicate  that  for  sufficiently  large  m,  the  convergence 
of  p",P2,  ...,Pp  to  Pi,P2,  "mPp  is  virtually  instantaneous,  and  that  (115)  is  a  nearly  perfect 
approximation  to  pg.  The  numerical  values  of  pi,P2*...,Pp  can  be  found  for  various  values  of  p 
in  Section  7.4.  Also,  note  that  cq  does  not  depend  on  p,  and  its  numerical  value  (to  16  digits) 
is  -.9189385332046727. 


The  proof  of  the  following  theorem  is  similar  to  the  proof  of  Theorem  4.6. 

Theorem  4.13  Suppose  that  5(ar)  =  |a:l^,  with  A  a  real  number  such  that  0  <  |A|  <  1,  and 
(110)  has  a  solution  p"  =  (Po5Pi»P2»'”>Pp'  Then  for  all  n  >  2p,  the  quadrature  weights 
PoiPiyp2i-'->Pp  independent  of  n. 


4.5  Corrected  trapezoidal  rules  for  other  singularities 

In  the  preceding  sections  quadrature  formulae  are  provided  for  singular  functions  of  the  form 

f{x)  =  d>(x)s{x)  + xp{x),  (116) 


and 

fix)  =  <Pix)s{x),  (117) 

where  the  singularity  s(x)  is  of  the  form  log{\x\),  or  x^  (0  <  |A|  <  1).  Obviously  the  procedure 
developed  in  the  preceding  sections  can  be  applied  to  other  singularities.  As  an  example,  we 
construct  a  quadrature  formula  to  approximate  the  definite  integral, 


where  /  is  of  the  form  (117), 


with  o  >  0,  and  ^(x)  €  c*[-a 


a(x) 


1 

■n/o*  —  X*’ 


kh,a  +  kh]  and  even  (i.e.,  ^(-x)  =  ^(a:)). 


(118) 


(119) 
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Remark  4.3  The  choice  of  the  singularity  (119)  is  dictated  by  the  frequency  with  which  it  is 
encountered  in  the  numerical  solution  of  partial  differential  equations,  in  signal  processing,  and 
other  areas.  Otherwise,  almost  any  integrable,  monotone  singularity  could  have  been  chosen. 

We  define  the  corrected  trapezoidal  rule  T”n  by  the  formula 

7’;n(/)=  J2  +  »'”/(!'.),  (120) 

i=-(n-2)  <=1 

where  h  =  a/{n  -  1),  Xj  =  jh,  yi  =  a  -  hi  for  1  <  i  <  k/2,  and  yi  =  a  +  h{i  —  k/2)  for 
k/2  +  1  <  i  <  k.  We  will  use  the  expression  T^n  with  appropriately  chosen  i/”  as  quadrature 
formulae  for  functions  of  the  form  (117),  and  the  following  construction  provides  a  tool  for 
finding  j/",  so  that  the  rule  is  of  order  2k  -  2,  i.e.,  there  exists  a  real  c  >  0  such  that 

I  ?’;«(/)-  r  fix)dx\<^.  (121) 

J —a  ^ 


For  an  even  integer  k  >  2,  we  will  consider  the  following  system  of  linear  algebraic  equations 
with  respect  to  the  unknowns  1/",  with  j  =  1,2,  ...,fc: 


•/-»  Va^-x^  i=-(n-2)  -  xJ 


(122) 


with  h  =  Xj  =  jh,  yj  =  a  -  hj  for  all  1  <  j  <  k/2,  and  j/j  =  a  +  h{j  -  k/2)  for  all 
^/2  +  1  <  J  <  A.  It  is  easy  to  see  that  the  linear  system  (122)  is  independent  of  the  length 
of  the  interval  a,  and  the  unknowns  "•>*'?  can  be  determined  by  solving  the  system  of 

equations 


*  /-i  Ur? 

"  L  VT^*‘  ■ 


). 


(123) 


with  h  =  Xj  =  jh,  yj  =  I  -  hj  for  all  1  <  j  <  k/2,  and  =  1  +  h{j  -  k/2)  for  all 
k/2  +  1  <  j  <  k. 

The  proof  of  the  following  theorem  is  quite  similar  to  the  proof  of  Theorem  4.1,  and  is  omitted. 


Theorem  4.14  Suppose  that  for  some  a  >  0,  f{x)  =  ^  ^  c*[-a  —  kh,a  +  kh]. 

Then  there  exists  such  c  >  0  that 


|r;.(/)-  l<;^ 


(124) 


for  all  sufficiently  large  n. 

The  authors  have  been  unable  to  construct  a  quadrature  rule  for  singularities  of  the  form 
(119),  which  is  independent  of  the  number  n  of  points  used  in  the  uncorrected  trapezoidal  rule. 
However,  this  is  a  relatively  minor  deficiency  since  the  weights  in  such  cases  can  be  precomputed 
and  stored. 
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5  Numerical  Results 


# 


9 


Algorithms  have  been  implemented  for  the  construction  of  the  quadratures  T^m, 

The  correction  coefficients  /3"*  are  calculated  using  (25),  and  (52).  In  the  tables  in  Section 
7.1  the  correction  coefficients  for  orders  of  convergence  upto  43  are  tabulated.  In  Table  1, 
convergence  results  are  presented  for  some  of  the  rules  T^m-  Column  1  of  this  table  contains 
the  number  of  nodes  discretizing  the  interval  [0, 1]  was  discretized.  In  column  2  are  the  relative 
errors  of  the  standard  1-sided  4th  order  corrected  trapezoidal  rule,  given  here  for  comparison. 
Columns  3-9  contain  the  relative  errors  for  the  rule  Tpm  for  various  orders  of  convergence  m. 
In  all  cases  the  integrand  was  of  the  form 

f(x)  =  sin(200a:)  +  cos(201x).  (125) 


The  quadrature  weights  for  the  rules  T\pmi  T'^kpmi  T^kpm,  and  r”n  are  all  obtained  as  solutions 
of  linear  systems,  and  it  is  easy  to  see  that  the  linear  systems  used  for  determining  these  weights 
(see,  for  example  (70),  (71))  are  very  ill-conditioned.  In  order  to  combat  the  high  condition 
number,  all  systems  were  solved  using  the  mathematical  package  MAPLE  using  200  significant 
digits. 

In  order  to  evaluate  the  coefficients  7  for  singularities  of  the  form  s{x)  =  |x|-^  or  s(x)  =  /o5(|i|), 
we  start  with  the  right-end  corrected  trapezoidal  rule  Tpm  of  order  40.  Under  these  conditions. 


I  l<  <126) 

71 

for  all  -fc  <  i  <  fcjJfc  3^  0  (see  Theorem  4.6)  and  for  reasonable  k,  the  convergence  of  aj*  to 
7i  is  almost  instantaneous.  The  construction  of  the  quadrature  weights  m  is  performed  in  a 
similar  manner.  In  Section  7.2  the  coefficients  7,  are  listed  for  the  singularities  lop(|x|),  |x|i, 
lx(“a,  |x|^,  1*1“^,  |®|~»  and  for  the  same  singularities,  the  quadrature  weights  m  are  listed 
in  Section  7.3.  In  Table  2,  convergence  results  are  presented  for  some  of  the  rules  TJ^kpm  for 
various  singularities.  Column  1  of  this  table  contains  the  number  of  nodes  in  the  discretization 
of  the  interval  [0, 1].  In  Table  3,  convergence  results  are  presented  for  some  of  quadrature  rules 
Tpkpm  for  various  singularities.  Column  1  of  this  table  contains  the  number  of  nodes  in  the 
discretization  of  the  interval  [—1, 1].  In  all  cases  the  integrand  was  of  the  form 


/(x)  =  (sin(20i)  4-  cos(21x))  -1-  (sin(23x)  -f-  cos(22x))s(x),  (127) 


and  the  order  of  convergence  used  was  10. 

Finally,  algorithms  have  been  implemented  for  evaluating  quadratures  to  integrate  func¬ 
tions  of  the  form 

fix)  =  <^(x)/off(|xl).  (128) 

The  quadrature  weights  are  obtained  by  solving  the  linear  system  (110).  Note  that  the  quadra¬ 
ture  weights  are  independent  of  the  discretization  h,  except  for  the  first  weight  po  which  is 
calculated  using  the  formula  (115).  Presented  in  Table  4  are  convergence  results  for  integrating 
functions  of  the  form  (128)  where, 

^(x)  =  sin(200x)  +  cos(201i).  (129) 
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Column  1  shows  the  number  of  nodes  in  the  discretization  of  the  interval  [-1,1].  Columns  3-6 
show  the  relative  errors  for  the  various  orders  of  convergence  m  as  shown. 


r 


Table  1:  Convergence  of  quadrature  rules  Tpm  for  non-singular  functions 


N 

k=4 

m=3 

in=9 

m=:15 

in=21 

in=27 

xn=33 

m=39 

20 

.230E-01 

.112E-01 

.131E-01 

.136E-01 

.138E-01 

.138E-01 

.138E-01 

.138E-01 

40 

.132E-01 

.120E-01 

.122E^-01 

.122E-01 

.122E-01 

.122E-01 

.122E-01 

.122E-01 

80 

.457E-02 

.108E-02 

.654E-03 

.430E-03 

.292E-03 

.202E-03 

.142E-03 

.lOOE-03 

160 

.216E-03 

.804E-04 

.223E-05 

.743E-07 

.264E-08 

.972E-10 

.365E-11 

.139E-12 

320 

.310E-05 

.522E-05 

.292E-08 

.199E-11 

.116E-14 

.304E-15 

.306E-15 

.306E-15 

640 

.191E-06 

.328E-06 

.304E-11 

.105E-15 

.703E-16 

.703E-16 

.703E-16 

.703E-16 

1280 

.239E-07 

.205E-07 

.266E-14 

.345E-15 

.346E-15 

.346E-15 

.346E-15 

.346E-15 

Table  2:  Convergence  of  quadrature  rules  for  singular  functions  (10th  order) 


N 

Ml*l) 

|i| 

!ili 

IllT 

40 

0.29128E-03 

0.25056E-04 

0.11650E-02 

0.42S10E-04 

0.53715E-03 

80 

0.72599E-07 

0.30493E-07 

0.98819E-06 

0.53217E-07 

0.52449E-06 

160 

0.56928E-10 

0.17499E-10 

0.10903E-08 

0.32715E-10 

0.49582E-09 

320 

0.65586E-13 

0.59119E-14 

0.76827E-12 

0.12962E-13 

0.31491E-12 

640 

0.18596 E-14 

0.16376E-14 

0.66613E-15 

0.17208E-14 

0.13878E-14 

Table  3:  Convergence  of  quadrature  rules  T^kpm  for  singular  functions  (10th  order) 


N 

/09(|i|) 

l*l» 

\z\^ 

li|i 

IzlT 

0.57489E-03 

0.49592E-04 

0.23137E-02 

0.84150E-04 

0.10655E-02 

80 

0.14438E-06 

0.60500E-07 

0.19680E-05 

0.10563E-06 

0.10436E-05 

160 

0.11348E-09 

0.34867E-10 

0.21762E-08 

0.65197E-10 

0.98921  E-09 

320 

0.13357E-12 

0.13614E-13 

0.15360E-11 

0.28103E-13 

0.62927E-12 

640 

0.61062E-15 

0.16237E-14 

0.42188E-14 

0.16237E-14 

0.50515E-14 

Table  4:  Convergence  of  the  quadrature  rule  T^kpm  for  functions  /(i)  =  <t>{x)log{\x\) 


N 

m=:3 

in=9 

in=lS 

in=21 

m=27 

in=33 

II 

B 

40 

.546E-01 

.S36E-01 

.536E-01 

.536E-01 

.536E-01 

.536E-01 

.536E-01 

80 

.291E-03 

.764E-03 

.265E-03 

.129E-03 

.640E-04 

.300E-04 

.107E-04 

160 

.282E-03 

.241E-04 

.209E-05 

.255E-08 

.482E-09 

.125E-11 

.143E-13 

320 

.437E-04 

.190E-04 

.912E-06 

.392E-09 

.162E-09 

.294E-12 

.147E-14 

640 

.573E-05 

.315E-05 

.468E-06 

.166E.07 

.108E-09 

.583E-13 

.119E-14 
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6  Generalizations  and  Conclusions 

A  group  of  algorithms  has  been  presented  for  the  construction  of  high-order  corrected  trape¬ 
zoidal  rules  for  functions  with  various  types  of  singularities,  both  end-point  and  in  the  middle  of 
the  interval  of  integration.  In  many  cases,  the  corrected  rule  can  have  effectively  an  arbitrarily 
high  order,  without  the  attendent  growths  of  correction  weights.  The  drawback  of  the  approach 
is  the  need  for  the  integrand  to  be  available  in  a  small  area  outside  the  interval  of  integration, 
whenever  the  singularity  being  corrected  is  on  one  of  the  ends  of  that  interval. 

The  algorithm  of  the  present  paper  admits  several  straightforward  generalizations. 


1.  There  are  classes  of  singularities  not  covered  by  this  paper  for  which  some  versions  of 
Theorem  4.1  can  be  fairly  easily  proven. 

2.  The  quadratures  can  be  easily  modified  to  handle  functions  of  the  form 


m 

fix)  =  rp(x)  +  Mx)  ■  »i(x), 

i=l 


(130) 


where  axe  smooth  functions,  and  si,  s^, ...,  Sm  are  several  different  singularities. 

3.  Quadrature  rules  developed  of  this  paper  have  fairly  obvious  analogues  in  two  and  three 
dimensions.  However,  the  proofs  of  the  multidimensional  versions  of  the  theorems  in  this  paper 
are  somewhat  more  involved  than  those  of  their  one  dimesional  counterparts.  These  results  will 
be  reported  at  a  later  date. 

4.  High-order  corrected  trapezoidal  rules  can  be  used  to  approximate  integrals 


by  rewriting  the  integral  as 


f  cosia  cos  6)d6 
Jo 

r  cosjx) 


(131) 

(132) 


and  using  the  quadrature  rule  TVn  defined  in  (120).  This  rule  proves  to  be  of  fundamental 
importance  in  the  development  of  the  Fast  Bessel  Ttansform  (see,  for  example  [10]). 

5.  Integral  equations  of  the  form 


ft 

I  a{w)logi\z  —  w\)dsru  =  C 
Jo 


(133) 


are  encountered  in  the  study  of  partial  differential  equations  (see,  for  example  [11]).  In  order  to 
apply  the  Nystrom  algorithm  to  the  integral  equation  (133),  the  left-hand  side  is  decomposed 
into  a  sum 


/  ffiw)logi\z  -  u;|)ds,j,  =  /(r)  +  J(z), 
Jo 

where  the  integral  operators  I  and  J  are  defined  by  the  formulae 

tL 


Hz)  =  I  ciw)logi\  \)ds^. 


(134) 


(135) 
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J{z)=  I  a{w)log{\‘y  (136) 

Now,  the  integral  operator  I  can  be  discretized  by  the  uncorrected  trapezoidal  rule  and  the 
operator  J  can  be  discretized  by  the  corrected  trapezoidal  rule  TJ),  defined  in  (109)  to  a  rapidly 
convergent  finite-dimesional  approximation  to  (133). 
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7  Correction  weights  for  Non-singular  and  Singular  functions 

% 

7.1  Quadrature  Weights  for  Non-singular  Functions 


/: 


f{T)dx 


7?-(/) 


=  r„(/)  +  fc  ^  {fib  +  kh)  -  f(h  -  kh)  -  f{a  +  kh)  +  /(a  -  kh))fiT. 

Jksl 


wmmmmBEEk 

m  =  S 

m  =  7 

0.4166666666666667D-01 

0.5694444444444444D-01 

-0.7638888888888889D-02 

0.648396164021 1640D-01 
-0.1395502645502646D-01 
0.1579034391534392D-02 

m  =  11 

3 

II 

0.696S636022927690D-01 
-0.1877177028218695D-01 
0.3643353174603175D-02 
•0.34405313051 14638D-03 

0.7289995064734647D-01 

-0.2247873075998076D-01 

0.57285184433621930-02 

-0.96187987681043240-03 

0.7722834328737J06D-04 

0.75232409136737010-01 

-0.25394303871718930-01 

0.76722338511876380-02 

-0.1739366039940610O-02 

0.25392974399877510-03 

-0.17670140071140400-04 

m  =  15 

m  t=  17 

m  =  19 

0.76990174607492560-01 
-0.27737991166059670-01 
0.94299993219431970-02 
-0.25916159651554270-02 
0.52025784562840550-03 
-0.66838404987379850-04 
0.409735540968662 1 0-05 

0.78362263347846430-01 

-0.29658915402555080-01 

0.11001664606348530-01 

-0.34647633453806100-02 

0.85608376109962970-03 

-0.15319364039426610-03 

0.17530392028535590-04 

-0.95950261563206930-06 

0.79463018590824320-01 
-0.31 260013937795620-01 
0.12402625824684000-01 
-0.43268933258947500-02 
0.12409632166862990-02 
-0.27635506618200040-03 
0.44471953919602460-04 
-0.4581 897491 74 1 90 1 0-05 
0.22639967975686450-06 

m  =  21 

m  =  25 

0.80365661345810830-01 

-0.32613978070275400-01 

0. 1 3652438870049960-01 
-0.51601020228053840-02 
0.16575675651416160-02 
-0.43258169685274430-03 
0.87357695672355700-04 
-0. 12750610206552040-04 
0.1193747238089644D-05 
-0.5374 1 53 101 848776D-07 

0.81119247515189910-01 
-0.33773341407781680-01 
0.14770396374073870-01 
-0.59550940256668330-02 
0.20923288167064710-02 
-0.6167158739860944O-03 
0.14703080863223770-03 
-0.2710805091 8704100-04 
0.36165653582653040-05 
-0.31 01 2440087834590-06 
0.12819143492992910-07 

0.81757875072513670-01 
-0.34776898997861870-01 
0.15773953964154060-01 
-0.6707762218226974O-02 
0.25350748123300830-02 
-0.82333067194378020-03 
0.22315204998506930-03 
-0.48856977019513130-04 
0.82770495227243840-05 
-0.101 62583651 903280-05 
0.803623922532694 1 0-07 
-0.3070147670921659O-08 
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27 


m  =  39 


0.8444980879146044D-01 
-0.3922700061890783D-01 
0.2063059004248263D-01 
-0. 1085151806728575D-01 
0.5447289134690455D-02 
-0.2547701211583463D-02 
0.1093355313271473D-02 
-0.42557271 1931 7083D-03 
0.1487041779510615D-03 
-0.4617000028477128D-04 
0.1 25959581 0865357D-04 
-0.2980436293579194D-05 
0.6019365929090900D-06 
-0.1016414607443390D-06 
0.1395254130438025D-07 
-0.1495069020432704D-08 
0.11 72703286200068D-09 
-0.5987202298631028D-1 1 
0.1492744845851982D-12 


m  =  41 


0.84688618781915600-01 
-0.3963949060242128D-01 
0.21 1 1481 7414433200-01 
-0.11308843918572410-01 
0.58131498157197770-02 
-0.28009893753729940-02 
0. 1 2465790 1 72923000-02 
-0.50687508549377960-03 
0.186551 83460926720-03 
-0.6158941596033656D-04 
0.18067363670950930-04 
-0.46591630001931570-05 
0.1042814313838009D-05 
-0. 1 9939262963808130-06 
0.3 190683763 1 802300-07 
-0.41 549647726433780-08 
0.42279889475231400-09 
-0.31526741882446180-10 
0.15317566842788960-11 
-0.36381110518255210-13 


m  =  43 


0.84905823450735190-01 

-0.40017237852542320-01 

0.21563392273951940-01 

-0.11739475783710390-01 

0.6165108551649863O-02 

-0.30512711431454990-02 

0.14030051221501160-02 

-0.59317914334636790-03 

0.22862506281240400-03 

-0.79685428090717950-04 

0.24909918257751390-04 

-0.69211645164908300-05 

0.16914765133645480-05 

-0.3590633249061523D-06 

0.65171565812650440-07 

-0.99088637012225160-08 

0.1227209106806963O-08 

-0.11888350699245330-09 

0.84475005888211280-11 

-0.39148991177844680-12 

0.8877720031504791O-14 


7.2  Quadrature  Weights  7*  for  Singular  Functions 


rf(x)dx  « 

Jo 


0.75188123386400250+00 

-0.6032109664493744D+00 

0.10738668308721570+01 

-0.72253709828678500+00 


0.49111698029675020+00 

-0.3176980828356269O+00 

0.71410805711892340+00 

-0.3875269545800468D+00 


0.16351359417233530+01 

-0.15331151513609710+01 

0.21437194469404900+01 

-0.17457402373028730+01 


0.14201135710357900+01 

-0.31252877971788190+01 

0.2592853861401367D+01 

-0.76486987895843140+00 

0.20277260836205720+01 

-0.37302381487966240+01 

0.29141056431S0046D+01 

-0.83440333427390050+00 


0.89518545428760170+00 

-0.1631355661694529O+01 

0.12165280228991150+01 

-0.33189682911689870+00 

0.13232780978696490+01 

-0.19969978433419440+01 

0.13925132311121590+01 

-0.36725447201515240+00 


0.31924164003655870+01 
-0.83495190059975070+01 
0.76531 189087438080+01 
-0.241 572142601 38580+01 
0.41277319448148460+01 
-0.94315385700363980+01 
0.828551 9053356245  D+01 
-0.25620073052327220+01 


28 


29 


0.553409i724301567D+00 

-0.3866961728429464D-i-00 

0.8032238407479816D+00 

-0.4699368403351921D+00 


=  2 


0.1181425202719417D+01 
-0.1060178333186577D+01 
0. 1 61 31 04391 254726D+01 
-0.1234351260787565D+01 


0.9469239981678674D+01 

-0.9440762908621185D+01 

0.1027199835611538D+02 

-0.9800475429172870D+01 


1 _ _ 

-1 

-2 

-3 

-4 

1 

2 

3 

4 

0.102083207138862504-01 

-0.198318610254488504-01 

0.153338124322483104-01 

-0.429839218127034704-00 

0.149833181703408204-01 

-0.241269929364636904-01 

0.174707110362580304-01 

-0.47389162095505300400 

0.22824861998852230401 
-0.56508138767703680401 
0.50151764926778740401 
-0.15496987012608240  401 
0.30836433412134590  401 
-0.65323932485360340401 
0.55143295626359260401 
-0.16627297698452560401 

0.18872991401279020402 
-0.55335856573322430402 
0.54466695681 138020402 
-0.1 7982569314390280402 
0.20317027997461520402 
-0.57223484572975360402 
0.55656419000928460402 
-0.18271223620118950402 

0.1453673785846622D+01 
-0.4645097217879781D+01 
0.713468108S431341D+01 
•0.S928340311S44518D+01 
0.2571788299099303D+01 
.0.4588965281604129D+00 
0.2135823382632839D-f01 
-0.5636852769577275D+01 
0.81094431 12743051D+01 
-0.6522597930471829D+01 
0.2775248991033700D+01 
-0.4888738991530396D+00 


0.3345150279872830D+01 
-0.1359274618S99862D+02 
0.2395745376553084D+02 
-0.2193610831631847D+02 
0.1025817941642386D+02 
•0.1946039989983884D+01 
0.4476264293482232D+01 
-0.1S61643865091390D+02 
0.2622632514046215D-I-02 
-0.2345727234258157D+02 
0. 1 081 90921 681 1830D-I-02 
-0.2033859578093758D+01 


0.2823540877500425D+02 

-0.1374572775395727D+03 

0.2701409207262959D+03 

-0.2668390934875709D+03 

0.1321842198719930D+03 

-0.2624585302793210D+02 

0.3025125337714398D+02 

-0.1418090842954654D+03 

0.2756109143783281D+03 

-0.2708086499560764D+03 

0.1337381187138368D+03 

-0.2650087753598454D+02 


s(i)  =  X 


0.1866196808675184D-(-01 
-0.8307135206229368D+01 
0.1909144688191794D+02 
-0.2636153756127239D+02 
0.2279974850816623D+02 
-0.1215894117169713D+02 
0.3671126621978929D+01 
-0.4816958588438264D+00 
0.2736714477854559D+01 
-0.1 004886340865 1 74D-I-02 
0.2164374286136325D-I-02 
-0.2894336017656429D-1-02 
0. 245606865484754404-02 
-0. 1293400965739323D-I-02 
0.3870327870200545D-I-01 
-0.S044475379801205D-f-00 


a(x)  =  x~ 


k=  8 


0.4383819645513359D-f01 
-0.2478753951 1 12947D-J-02 
0.6535630997043997D4-02 
-0.9943323854718145D+02 
0. 9269762740745608D+02 
-0.52554488204227320-1-02 
0.16709325789196860-1-02 
-O.22927745642297460-«-Ol 
0.581 91 435971 649600-f  01 
-0.283373291149343104-02 
0.713000686392113204-02 
-0.106050880126925604-03 
0,975610176847412404-02 
-0.5482984235409941 04-02 
0.173251010867274604-02 
-0.236632139772391104-01 


0.3756991 225931 81 304-02 
-0.255656820049079804-03 
0.753156064006868204-03 
-0.123906065368688704-04 
0. 1226707735965091 04-04 
-0.730098332404368404-03 
0.241736444448029404-03 
-0.343377107423119104-02 
0.401 1 57503348349504-02 
-0.263313343863456904-03 
0.767575136075255204-03 
-0.125648702569937504-04 
0.124034129581003104-04 
-0.736805543681391004-03 
0.243629021453027904-03 
-0.34571 9302255855304-02 
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k=  10 


0.2264788460960479D+01 

0.54054546335160520-1-01 

0.46883768289745560-1-02 

-0. 1292939169279749D4-02 

-0.3917131575943378O-f02 

-0.40983301861233700-1-03 

0.3957101757244672D4-02 

0.13752207611158520-1-03 

0.1609252630383255O-f04 

-0.7639785056816069D4-02 

-0.29252186647286950-1-03 

-0.37053009775985810-1-04 

-5 

0.98982375845036270-1-02 

0.40850243893952150-1-03 

0.55008520533332450-1-04 

-6 

-0.8785457780822613D-f02 

-0.38544638637505400-1-03 

-0.54548359992901400-1-04 

-7 

0.5297228202020211D-I-02 

0.24472818048521790-1-03 

0.36110062477032520-4-04 

-8 

-0.2081786990425939D-1-02 

-0.10057549373541660-1-03 

-0.15382377271423390-1-04 

-9 

0.4823072180507942D-1-01 

0.24238208038636850-1-02 

0.38253464966205660-1-03 

-10 

-0.50078745884467410-1-00 

-0.26069872827045750-1-01 

-0.42306114818517930+02 

1 

0.331 16208887872880-1-01 

0.71265780202799180-1-01 

0.49930633080662240+02 

-0.15591265293058690-1-02 

-0.44600419155155310-1-02 

-0.42157854458783690+03 

0.44752777942639230-1-02 

0.149613880576318604-03 

0.16387301739716230+04 

-0.83719548276891600-1-02 

-0.31133829920638690-1-03 

-0.37551629380686000+04 

5 

0.10645933091384000-1-03 

0.42921423110749820-1-03 

0.55593520581381000+04 

6 

-0.93330160619516560-1-02 

-0.40157257212274240-1-03 

-0.55027899293502870+04 

7 

0.55782099287843460-1-02 

0.2534431424358416O-f03 

0.36380609038511120+04 

8 

-0.21778940784503470-1-02 

-0.10369302010580000-1-03 

-0.15482791409953160+04 

9 

0.50201723812649580-f01 

0.24903333917958600-1-02 

0.38474701600188970+03 

10 

-0.51914508726977670-1-00 

-0.26711640508111780-1-01 

-0.42525743950987800+02 

7.3  Quadrature  Weights  /i*  for  Singular  Functions 


=  (/)  +  Tl/J"  (/)  +  #*;(/(*;)  + /(*->)) 


H! 

a(x)  =  log(x) 

s(z)  =  x> 

1 

2 

0.18257480647361590+01 
-0. 13257480647361590+01 

0.12052250374156740+01 

-0.70522503741567370+00 

0.37788553886638430+01 

-0.32788553886638430+01 

II 

1 

0.34478396546563620+01 

-0.68555259459754430+01 

0.55069595045514130+01 

-0.15992732132323320+01 

0.2218463552157251O+01 

-0.3628353505036473O+01 

0.26090412540112730+01 

-0.69915130113205120+00 

0.73201483451804340+01 

-0.1778105757603391O+02 

0.15938637962100050+02 

-0.49777287312465800+01 

II 

1 

0.49673629782877580+01 
-0.16205015048591260+02 
0.25851537618326390+02 
-0. 2222599466791883 0+02 
0.99301049980375390+01 
-0.18179958781415940+01 

0.31437306974379690+01 

-0.84518976058583290+01 

0.12083574361193640+02 

-0.96175706597041530+01 

0.40477322199084100+01 

-0.70556901297753240+00 

0.10736553147151140+02 

-0.42999806712146430+02 

0.7669950880443822O+02 

-0.71150124092010390+02 

0.33666381705342880+02 

-0.64525128527754170+01 
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0.64226100151584150+01 

-0.2917804569054941D+02 

0.6966653255738346O+02 

-0.9953245609468264D+02 

0.S864202397339461D+02 

-0.4848666190307500D+02 

0.14969634046024750+02 

-0.20036369036541830+01 


0.40151621208378910+01 

-0.15031948829934550+02 

0.32174891825488320+02 

-0.42584652819626590+02 

0.35755327013269900+02 

-0.18640786941749750+02 

0.55260187983314370+01 

-0.7140111666166550D+00 


0.14073427412263150+02 

-0.78654272014335400+02 

0.2101259829107628D+03 

-0.32416246252361000+03 

0.30600514496632690+03 

-0.17540699865587930+03 

0.56312989573782880+02 

-0.77938116693110270+01 


0.78324320205687790+01 

-0.45651616703747490+02 

0.14521688463546780+03 

-0.29013483028863790+03 

0.38708621625799000+03 

-0.35238213835706800+03 

0.21724215475193420+03 

-0.87077960873829890+02 

0.20535842660726350+02 

-0.21669841034038230+01 


=  5 


0.48489969696179590+01 
-0.23281314990198370+02 
0.66396083214006050+02 
-0.12292552699668700+03 
0 . 1 5467087356924 190+03 
-0.13425144850091400+03 
0.79503700214730130+02 
-0.30777567756058370+02 
0.70387973129601030+01 
-0.72259303669838690+00 


0.1735350975431 1530+02 
-0.12456287138310560+03 
0.44346408794764410+03 
-0.95699764679509340+03 
0.13537838609265590+04 
-0. 12919502347499440+04 
0.82855034879631690+03 
-0.34356260419817710+03 
0.83464206930537340+02 
-0.9042657229049132D+01 


k=  1 


0.13566330131781380+01 

-0.85663301317813840+00 


0.27945295939741420+01 

-0.22945295939741420+01 


0.25191638884227070+01 

-0.43958853961912530+01 

0.32804523468506340+01 

-0.90373083908208770+00 

0.53661295410986820+01 

-0.1218320712530640O+02 

0.10529506055313800+02 

-0.32124284711060800+01 

k=  i 

0.35894971 684794600+01 
-0.10281949987457060+02 
0.15244124198174390+02 
-0.12450938242016350+02 
0.53470372901330030+01 
-0.94777042731345250+00 

0.78214145733550620+01 

-0.29209184836912520+02 

0.50183778905992990+02 

-0.45393380658900040+02 

0.21077271584542160+02 

-0.39798995680776420+01 

A 

0.460291 1 2865297440+01 
-0.18355998614881 100+02 
0.40735189743281 190+02 
-0.55304897737836680+02 
0.4736043505664168D+02 
-0.25092950829090350+02 
0.75414544921794740+01 
-0.98614339682394690+00 

0.10202963242678320+02 

-0.53124868626063780+02 

0.13665637860965130+03 

-0.20548411867410710+03 

0.19025864509219730+03 

-0.10738433055832670+03 

0.34034426875924320+02 

-0.46590959619536570+01 
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It 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

0.5576409349747767D+01 
-0.2852065698585619D+02 
0.8432379551508595D+02 
-0. 1 601 173988450523D+03 
0.2054417067588763D+03 
-0.1811847384277427D+03 

0. 1 08754381 3080456D+03 
-0.4259681068876286D+02 
0.9843244561772901D+01 
-0.1019932546114451D+01 

0.12S3203265379597D+02 

-0.8377173491458909D+02 

0.2871359566879038D+03 

-0.6038601656792564D+03 

0.8377166700470196D+03 

-0.7870189584977965D+03 

0.4981713229210595D+03 

-0.2042685138412166D+03 

0.4914154195659545D+02 

-0.5278151333515754D+01 

7.4  Quadrature  Weights  for  Singular  Functions 


Tr4.»  (/)  +  TI;,-.  {/)  +  A  53  PAH^j)  + 

J=0 


Noie:  po  is  given  for  A  =  0.01,  For  any  other  h,  the  following  formula  is  used  to  calculate  po. 


p 

po  =  (-.9189385332046727417803  +  0.5log(h))  -  53  Pj 


(137) 


m  =  3 

m  =  5 

m  =  7 

-0.3221523626198730D+01 

-0.3191075169140337D-(-01 

-0.3044845705839327D-01 

-0.318146710201317lD-(-01 

-0.4325921322794724D-01 

0.3202689042388491D-02 

m  =  9 

m  =  11 

m  =  13 

-0.3176811195217937D+01 
-0.5024307342079858 D-01 
0.5996233119529027D-02 
-0.4655906795234226D-03 

-0.3174071153542312D-I-01 
-0.5462714010179898D-01 
0.81 88266460029230D-02 
•0.1091885919666338D-02 
0.7828690501786438D-04 

-0.3172268092036274D-(-01 

-0.5763224261186158D-01 

0.990546789435U714D-02 

-0.1735836457536894D-02 

0.2213870245446548D-03 

-0.1431001195267904D-04 

m  =  15 

m=  17 

m  =  19 

-0.3170992165916170D-f01 
-0.5981954453204053D-01 
0.1127253159446256 D-01 
-0.2343420324253269D-02 
0.403662184559S672D-03 
-0.4745095013720858D-04 
0.2761744848710794D-05 

-0.3170041916020681D-I-01 
-0.6148248184914681D-01 
0.1238115647253341D-01 
-0. 2897732763288696D-02 
0.6052303442088134D-03 
-0.97842990049S2013D-04 
0.1051436637368180D-04 
-0.5537586803550720D-06 

-0.3169306861514305D-I-01 
-0.6278924541603596D-01 
0.1329589096935583D-01 
-0.3396678852464559D-02 
0.8131245480320894D-03 
-0.1618104373797589D-03 
0.2422167651587582D-04 
-0.2381400032647607D-05 
0.114227584S182835D-06  | 
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m  =  21 


m  =  25 


-0.3168721384920306D+0] 
-0.6384310328523506D-01 
0. 1406233305604607D-01 
-0.3843770069700536D-02 
0.1019474340602541D-02 
-0. 235506791 8692057D-03 
0.4387403771306165D-04 
-0.60662177571 19951D-05 
0.5477355521032650D-06 
-0.2408377597694342D-07 


-0.3168244070581230D+01 
-0.6471094753809971D-01 
0.1471321624569456 D-01 
-0.4244313571022683 D-02 
0.1219746091263614D-02 
-0.3156154921336350D-03 
0.6890800654569580D-04 
-0. 1 19S656336479857D-04 
0.1S29459820049702D-05 
-0.1274231726028842D-06 
0.5166969831297037D-08 


-0.3167847493678468D-I-01 
-0.6543800519316396D-01 
0.1527249136497475D-01 
-0.4603847576274231D-02 
0.1411497560731106D-02 
-0.3995067600256630D-03 
0.985166893311174SD-04 
-0.2018119747186014D-04 
0.3260961737325821D-05 
-0.3871484601943021D-06 
0.2990271 150667017D-07 
-0.1124351894335142D-08 


_ m  =  27 

-0.3167512774719844D-J-01 
-0.6605594788600917D-01 
0.1575801776649599D-01 
-0.4927S31843955059D-02 
0.1593569961301572D-02 
-0.48S1878897058822D-03 
0.1318371286512027D-03 
-0.3070344146767653D-04 
0.S891S22736279919D-05 
-0.88820769809032060-06 
0.9822897121976360D-07 
-0.7065765782430224D-08 
0. 24755891 2003961 70-09 


_ m  =  29 

-0.31672264928891640-1-01 
-0.66587614142985770-01 
0.16183350772077270-01 
-0.52199482852921880-02 
0.17655796326763540-02 
-0.57119272539327300-03 
0. 1 68049691 04589350-03 
-0.43377838305818330-04 
0.95127789757490050-05 
-0.17112204797878400-05 
0.24136162890628880-06 
-0.24957347993245870-07 
0. 1 6788854888692140-08 
-0.55051022187125070-10 


_ m  =  31 

-0.31669788467725300-1-01 

-0.67049886894035770-01 

0.16558947382305380-01 

-0.54850753042767410-02 

0.19276016998335810-02 

-0.65646749758128730-03 

0.20642333853049990-03 

-0.57996370680906480-04 

0.14164130186004330-04 

-0.29246164476805320-05 

0.49415245555059960-06 

-0.65403880256335600-07 

0.63457930576872600-08 

-0.40074787913661000-09 

0.12346316319624460-10 


m  =  33 

-0.31667625116197080-1-01 

-0.67455515305576880-01 

0.16892994309456880-01 

-0.57263314183306020-02 

0.20799739823939140-02 

-0.74027225298947030-03 

0.24633036491534910-03 

-0.74321972383799320-04 

0. 19842600343532270-04 
-0.45808369102928480-05 
0.8«6453665775555O-06 
-0.14184482468459630-06 
0.1767037741742959O-07 
-0.16140962033947170-08 
0.9602551109604562O-10 
-0.27893064925473720-11 

_ L _ 

-0.31665719037402870-1-01 

-0.67814306608015600-01 

0.17191987061489160-01 

-0.59466418671964910-02 

0.22231757741567420-02 

-0.82210184828251480-03 

0.28724516256187130-03 

-0.92111014838808970-04 

0.26513491264160890-04 

-0.67155220048940090-05 

0. 14663682766624830-05 
-0.26956102692569140-06 
0.40476842103339430-07 
-0.47598154704167640-08 
0.41059743779825030-09 
-0.23084269505592830-10 
0.6342175941576707O-12 

-0.31664026923256750-1-01 
-0.68133928168950600-01 
0.17461142060171270-01 
-0.61485081 162080720-02 
0.2357753273497796O-02 
-0.90162491607495610-03 
0.32873545880140590-03 
-0.11112740061526230-03 
0.34120045574742230-04 
-0.93485600354798580-05 
0.22465276931323640-05 
-0.464600881 043161 70-06 
0.80829915369022930-07 
-0.1 148532768136401 0-07 
0. 1 2784054650] 72500-08 
-0.10444127205737410-09 
0.55649450215383520-11 
-0.14502139492296120-12 

m  =  39 

II 

e 

B 

tl 

-0.3I6625146677508lD-h01 
-0.68420460791 12800D-0I 
0.I770469478902206D-01 
-0.6334072100094389D-02 
0.2484274171602103D-02 
-0.9786376366601 862D-03 
0.3704506824517389D-03 
-0.1311 507079674222D-03 
0.42S9144483911754D-04 
-0.1248611531858182D-04 
0.3255027605557997D-05 
-0.7428077534364396D-06 
0.1457448522607878D-06 
-0.2404950901525398D-07 
0.3241558798437558D-08 
-0.3423992518658962D-09 
0.2656123735758442D-10 
-0.1344809S28411308D-11 
0.3332744815245408D-13 

-0.3166115504658825D+01 
-0.6867878881201404D-01 
0.1792611880692437D-01 
-0. 65051 72477564357D-02 
0.2603300521 146428D-02 
-0.1053029105125390D-02 
0.4121099047922528D-03 
-0. 1 51 98031 91 376791D-03 
0.5184904980367620 D-04 
-0.161230315546S844D-04 
0.45091366524809680-05 
-0. 1 1 1 90404675133310-05 
0.24283716557095330-06 
-0.45288452551852680-07 
0.7103184os-»l  309580-08 
-0.91028544268-104330-09 
0.91462516303229810-10 
-0.67532494409650900-11 
0.32567555153373960-12 
-0.76933711416127780-14 

-0.31661 155046588250-1-01 
-0.68678788812014040-01 
0.17926118806924370-01 
-0.65051 724775643570-02 
0.260330052 1 1 464280-02 
-0.10530291051253900-02 
0.41210990479225280-03 
-0.1 51 98031 91 376791 0-03 
0.51 849049803676200-04 
-0.16123031554658440-04 
0.45091366524809680-05 
-0.11190404675133310-05 
0.24283716557095330-06 
-0.45288452551852680-07 
0.7103184896000958O-08 
-0.91028544268404330-09 
0.91462516308229810-10 
-0.67532494409650900-11 
0.32567555153373960-12 
-0.76933711416127780-14 
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Errata: 


+  /(^l)  +  fi-Cn-l)  - 


1.  Equation  (1)  should  be  changed  to 

f  n.)d.  =  M  £  /(x<  j  +  . 

«=i  ^  24 

2.  Equation  (2)  should  be  changed  to 

^(/)  =  .(g  nx,) . 

J=I 

3.  The  weights  tabulated  in  Section  7  1  are  —3^  anA  ■jk  l  ,  , 

and  not  as  they  should  be. 


fi 


